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A group of non-uniform quantum lattice Hamiltonians in one dimension is introduced, which 
is related to the hyperbolic 1 + 1-dimensional space. The Hamiltonians contain only nearest 
neighbor interactions whose strength is proportional to cosh j A, where j is the lattice index and 
where A > is a deformation parameter. In the limit A ^ the Hamiltonians become uniform. 
Spacial translation of the deformed Hamiltonians is induced by the corner Hamiltonians. As 
a simple example, we investigate the ground state of the deformed 5 = 1/2 Heisenberg spin 
chain by use of the density matrix renormalization group (DMRG) method. It is shown that the 
ground state is dimerized when A is finite. Spin correlation function show exponential decay, 
and the boundary effect decreases with increasing A. 
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1. Introduction 

The density matrix renormalization group (DMRG) 
method has been applied to various problems in low di- 
mensional correlated physics. ^^"^ The method contains 
the block spin transformation in its formulation, but the 
relation with Wilson's renormalization group^ (RG) is 
not clear, since the hierarchy in energy scale is missing 
in the uniform lattice Hamiltonians. Recently Okunishi 
considered a group of half- infinite non- uniform systems, 
whose lattice Hamiltonians 

N N 

H{\)^Y. A^ ^..+i = E ^"ho+i (1-1) 

are sum of neighboring interaction terms ^j^^, where j 
is the lattice index and where -I- 1 is the system size.® 
He applied numerical RG method to this system, where 
the application is a one-parameter deformation — the 
exponential deformation — to the real space RG scheme 
introduced by Xiang.^^^ It is shown that the energy scale 
introduced by the deformation parameter A = e"^ > 1 
regularize the distribution of low energy excitations, even 
for the models that are gapless when A = in the large 
N limit. 

It is possible to introduce Okunishi's deformation 
scheme to the DMRG method for the purpose of regular- 
izing low energy excitations. A primitive way is to joint 
the deformed half-infinite Hamiltonians in Eq. (1.1) at 
the origin 

N 
j=-N 

to construct the whole system of size 2N + 2. The inter- 
action strength increases with \j\ toward the boundary 
of the system, in contrast to the smooth boundary con- 
dition proposed by Vekic and White. ^ The construction 
of the above Hamiltonian is, however, rather artificial in 
the point that one may choose arbitrary increasing func- 
tion of \j\ instead of e'-'l^. In this article we introduce a 



natural candidate 

N N 

J=-N j=-N 
N 

= cosh(jA) /i^ , (1.3) 

j=-N 

which keeps some aspects of translational invariance even 
when A > 0. It is possible to find a geometric interpreta- 
tion of H'^°^^{X) as a time boost in the hyperbolic 1-1-1- 
dimensional space, as discussed at the end of this article. 

In the next section we shortly review the exponen- 
tial deformation on the lattice Hamiltonian. Eigenvalue 
distribution is considered in the large system size limit. 
In §3 we explain the details of the deformed Hamilto- 
nian H""'"'^ (A) in Eq. (1.3). It is shown that the lattice 
translation is related to the deformed corner Hamilto- 
nian. Recursive construction of H'^"^^ (A) is also consid- 
ered. Ground state property of the deformed system is 
analyzed in §4 in the case of the 5=1/2 Heisenberg 
spin chain. It is shown that the spin correlation func- 
tion of the ground state decays exponentially when A is 
finite, and the state is dimerized. The effect of deforma- 
tion on the quantum entropy is observed. Conclusions 
are summarized in the last section. Quantum-classical 
correspondence of the deformed Hamiltonians and sev- 
eral conjectures are discussed. 

2. Exponential Deformation 

Consider a group of ID quantum Hamiltonians 

N 

^= E hj+i (2-1) 

on the lattice of size 2N -\- 2, where j+i represents 
the interaction between neighboring sites labeled by j 
and J + 1. A typical example is the antiferromagnetic 
Heisenberg spin Hamiltonian 

N 
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where J > represents the interaction parameter. Be- 
fore considering the deformed Hamiltonian H'^°^^{X) in 
Eq. (1.3), let us observe effects of the exponential defor- 
mation in Eq. (1.1). For latter convenience we treat the 
system whose linear size is 2N + 2. The exponentially 
deformed Hamiltonian is then written as 



jy'''^P(A) = 



N 

E 



(2.3) 



3 = -N 

whore the deformation parameter A is real and positive.^" 
When A = the above Hamiltonian H'^^p{X) coincides 
with the uniform Hamiltonian H in Eq. (2.1). 

It is known that the factor A = controls the eigen- 
value structure.^' ^ In order to observe the fact briefly, 
let us consider the infinite system size limit N oo. To 
simplify the discussion we assume that the ground state 
energy Eq is zero, and all other eigenvalues are positive. 
This assumption can be satisfied by adding appropriate 
constant to each neighboring interaction hj .j_^_i}^ 

Consider a right shift operation S that moves the lat- 
tice sites by one to the right direction. It is obvious that 
S^, the conjugate of S, represents the left shift operation, 
and therefore SS^ = S''''S' = 1 is satisfied. If we apply S 
to H'^^P{X) when the system size is infinite, we obtain 
the following relation 



5iI<'"P(A)5t = J2 '^'^ i^hj+i^^) 



(2.4) 



= £ e(^-i)^/i,.,.+i = e-^if-P(A). 



J = -OD 

As a result of translation the deformation parameter e^^ 
is modified to e^-'"^^^, and this modification can simply 
be expressed by multiplying the factor to H°^'^{X). 
This translation property in H^^p{X) restricts the eigen- 
value structure, which is obtained from the eigenvalue 
relation 



(2.5) 



If there is an eigenstate \'^) the shifted state S\^) is also 
an eigenstate, since we have the relation 

[SH^''P{X)S^S\^) = SH''''P{X)\-^) =ES\^), (2.6) 

and using the relation in Eq. (2.4) we can verify that 

[e-^ iJ'="P(A)] S\-^)^ES I*) (2.7) 

is satisfied. Thus if the eigenvalue E in Eq. (2.5) is posi- 
tive, there is a family of eigenvalues 

. . . , E, e-^ E, E, E, e^^ E 



? • • • ; 



(2.8) 



that are equidistant in logarithmic scale. Such a positive 
energy eigenstate \^) is not translationally invariant, and 
the orthogonality 



(^IS'I*) =0 



(2.9) 



It should be noted that presence of periodic eigenstates 
are not excluded. For example, if there is unique zero- 
energy eigenstate |$), it is translationally invariant. This 
is because Eq. (2.7) shows that S \^) is also the zero 
energy state. Thus we can say that if the zero-energy 
state is unique, it satisfies the translational invariance 



5|$) = |$). 



(2.10) 



As an extension one can consider digenerated case, where 
there are two zero-energy eigenstates and |$(,) that 
satisfies 



(2.11) 



This is the case when there is dimerization in the ground 
state. This degeneracy would be lifted by the effect of 
boundary when the system size 27V -I- 2 is finite. It is 
straightforward to extend the argument of degeneracy to 
trimerized state, etc. 

It is possible to consider various generalizations of 
H'^^P{X). As an example one can consider the deformed 
tight-binding Hamiltonian 

oo 



(2.12) 



for spinless lattice Fermions, where t represents the hop- 
ping parameter and where A the band gap. Since this 
Hamiltonian contains oscillating potential, the transla- 
tion period is 2-sitc when A = 0. Thus for this de- 
formed Hamiltonian H^^{X) one should modify the rela- 
tion Eq. (2.4) according to this period. It can be verified 
that all the one-particle states |^') satisfy the orthogo- 
nality in Eq. (2.9), and are represented by localized wave 
functions similar to wavelet basis function. The half-filled 
state 1$) has finite excitation gap, where 1$) is periodic 
and satisfies 5^ |$) = |$). When A = the one-particle 
eigenfunctions and energy spectrum is explained by the 
Bloch's theorem. It is not trivial how such an energy 
structure is destructed by the introduction of exponen- 
tial deformation. It is straightforward to generalize the 
exponential deformation to systems that contain inter- 
actions of longer range. 

3. Hyperbolic Deformation 

The eigenvalue distribution of H'^^p{X) explained in 
the last section prevents numerical study of the bulk 

property of the system around the center j = 0. This 
is because the energy scale in the left side of the system 
{j < 0) is smaller than that at the center, and to ap- 
ply the DMRG method to such system is difficult. This 
problem can be avoided if we take an average between 
fl'^^P(A) and H<'''p{-X) as 

if (A) = i [H'^'^P (A) + ""P ( - A)] 



JV 



is satisfied. 



^ coshjA /i^- . 



(3.1) 
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We call the deformation from H in Eq. (2.1) to i?'=°"'^(A) 
introduced here as the hyperbolic deformation in the fol- 
lowing. 

Let us extend the shift operation 5* and its conjugate 
S'^ to Hamiltonians of finite size systems. A natural way 
is to consider that the operation modifies the coefiicients 
of the neighboring interactions as follows 

cosh jA cosh(j - 1)A . (3.2) 

Then the shift operation on H''°^^{X) is defined as 

N 

SH^<'-^{\) = cosh(i - 1)A ft^. j.+i . (3.3) 

3=-N 

Taking the weighted difference between H'^^^^^X) and 
5^cosh(^-) S't we obtain the relation 

H'^osh|^^^ ^— S F™«'^(A) S't (3.4) 

cosh A 

N 

= tanh A ^ sinh jX hjj+i = tanh Ai?''"''(A) , 

j=-N 

where if''™^(A) introduced here represents deformed 
Hamiltonian of another type 

N 
j=-N 

which is decoupled at the origin j = 0. Similar to 

Eq. (3.4), the deformed Hamiltonian H™"'^(A) can be 
obtained from H^^'^^{X) by the following weighted differ- 
ence 

H^"^^\X) \— S H'''"'^(X) 5+ = tanh Ai/™"'^(A) . 

cosh A 

(3.6) 

The relations Eqs. (3.4) and (3.6) can be regarded as 
one parameter deformation to the translational invari- 
ance SHS^ = H, which is satisfied by the uniform Hamil- 
tonian in Eq. (2.1). 

Following the convention in the infinite system DMRG 
method, let us divide H'^°^^{X) into three parts 

H^-^\X) = H^iX) + /io,i + Hj,iX) , (3.7) 

where Hj^{X) and -ffjj(A) are defined as follows 

-1 

HlW = Yl coshjA /ij- 

j=-N 
N 

Hr{X) = 5^coshiA/i^.^.+i. (3.8) 

We also divide iJ^'"'^(A) in the same manner 

i/^inh(A) = CJA) + C«(A), (3.9) 

where Cl(A) and Cj^(A) are defined as follows 

-1 

j=-N 

N 

Cr{X) = ^sinhjA/i^-^-+i. (3.10) 
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These are deformations to the corner Hamiltonian,-^^' 
since in the limit A ^ we obtain the relation 

]im^^^=yjh...,. (3.11) 

A^o sinh A ^'^^^ ^ ' 

i=i 

We have shown the relation between H'^°^^{X) and 
^smh^^-j jtqj, ^j^g same system size 2iV + 2. We then focus 
on recursion relations, which connects systems of differ- 
ent sizes. Let us introduce Baxter's star notation^^' 

N 
N 

Q(A) = Esinh(j-l)A/i,-,.+i. (3.12) 
We then obtain recursion relation 

N 

Cr{X) = Y - + ^] (3.13) 

= cosh A C^(A) + sinh A [h^^^ + HrW] , 
and similarly we obtain 

N 

Hr{X) = Y ^°«^[(J' - + ^] ^-.i+i (3.14) 

= cosh A [/ii 2 + H^W] + sinh A Q(A) . 
If we introduce the double star notations 

N 

H*R*W = E^°^^(j-2)A^-,i+i 

N 

C*R*W = Esi^Mi-2)A/i,.,.+i, (3.15) 

we can decouple the recursion relations as follows 

iJ^(A) = cosh A /ii_2 - ^2,3 + 2 cosh A iJfl(A) - H^{X) 

CjiiX) = sinhA /ii_2 + 2coshA C^(A) -Q*(A). (3.16) 

These relations would be of use when one applies nu- 
merical renormalization group methods^' ^' ^' to the de- 
formed Hamiltonian H'^°^^[X) in order to obtain its eigen- 
states. 

4. Numerical Observations 

One might conjecture that the hyperbolic deformation 
violates uniform property of the system, since the bond 
interaction strength is modified. But for the ground state 
this intuition is not always true. For example, one can 
show that the valence bond solid (VBS) state of S' = 1 
spin chains is not violated by the hyperbolic (or even 
exponential) deformation. We observe another example, 
the ground state of the deformed 5=1/2 Heisenberg 
spin chain in this section. 

Figure 1 shows the nearest neighbor spin correlation 
function (s^ ^"j+i) calculated for the ground state of 400- 
site system when A = 0, 0.05, and 0.1. We keep m = 
130 states at most for the block spin variables in the 
calculation by the finite system DMRG method. When 
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Fig. 1. Nearest neiglibor spin correlation function (s? ^f+i) of the 
deformed 5=1/2 Heisenberg model. In all cases shown here the 
function contains even-odd oscillation, which decays very slowly 
only when A = 0. 



A = the correlation function show even-odd osciUation 
with respect to j, and the oscillation slowly decays from 
the boundary to the center of the system. It is known that 
the decay is in power low, which represents the gapless 
nature of the undeformed 5 = 1/2 Heisenberg chain. 
When A is finite, the oscillation is strongly stabilized, 
and the boundary effect disappears rapidly. 
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Fig. 2. Decay of the correlation function (— 1)^-'^^ (s^^ s^^j^) 
with respect to the distance 2j -|- 1. 



Figure 2 shows the correlation function | (s^^ •^f+i) I = 
(_X)2i+i ^gZ^ Sj+i) with respect to the distance 2j + 1. 
When A = 0.05 and A = 0.1 we observe exponential 
decay. The correlation length ^ obtained from the de- 
cay rate is almost inverse proportional to A as shown 
in Fig. 3, where ^A ~ 0.134 is satisfied. These calcu- 
lated results suggest that the hyperbolic deformation en- 
hances the local property of the system. To confirm this 
locality, we calculate the entanglement entropy. Figure 4 
shows the bipartite entropy 5 at the center of 400-site 
system. The value of S decrease exponentially with A, 
where S = 1.145 at the infinite A limit. 



Fig. 3. Correlation length ^ obtained from the spin correlation 
function in Fig. 2. 
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Fig. 4. Entanglement entropy 5 as a function of A. 



5. Conclusions and Discussions 

We have investigated the effect of hyperbolic defor- 
mation on ID quantum lattice Hamiltonians. Numeri- 
cal analysis on the deformed S* = 1/2 Heisenberg model 
shows that the deformation introduces dimerization in 
the ground state, and the local property is enhanced. 
The entangle entropy becomes finite even in the large 
system size limit when A > 0. 

Though the calculated system is one-dimensional it 
can possess a dimerized ground state, because shift of 
dimerized pattern introduces macroscopic increase of en- 
ergy expectation value; the dimer order might survive in 
finite temperature. In this sense the deformed system has 
property of higher dimensional systems. 

It would be interesting to consider whether the ground 
state is exactly represented by a matrix product state 
of finite matrix dimension in the infinite A limit. We 
conjecture that integer spin Heisenberg spin chains un- 
der strong hyperbolic deformation have such finite di- 
mensional matrix product ground states, if appropriate 
boundary conditions are imposed. 

The hyperbolic deformation can be used for scal- 
ing analysis of the ground state of undeformed system. 
The two parameter scaling proposed by Tagliacozzo et. 
g^j 14, 15 -^^j^gj-g |;jjg controllable parameters are the system 
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size N and the kept number of block spin states m, can 
be modified to the scahng analysis with respect to A and 
TO. More simply, if one keeps sufficient number of m in 
numerical calculation, the finite size scaling with respect 
to N can be replaced by the finite 1/A scaling. 

Let us cast our eye to the quantum-classical correspon- 
dence. If one considers the Trotter decomposition^^' of 
the deformed Hamiltonian H'=°''^{\), one finds that the 
Hamiltonian describes real or imaginary time boost in 
the hyperbolic 1 -|- 1-dimensional space, which has con- 
stant negative curvature. It is known that classical lattice 
models on the hyperbolic 2D space tend to show gaussian 
universality in their phase transition. This suggest 
that if H'^°^^{X) describes quantum phase transition of 
second order, it would subject to Gaussian universality 
class. Such a geometric interpretation may draw defor- 
mation of various type, such as spherical deformation 
iJ"'"(A). 

The authors thank to K. Okunishi for valuable dis- 
cussions and critical reading of this article. This work 
is partially supported a Grant-in-Aid for Scientific Re- 
search from Japanese Ministry of Education, Culture, 
Sports, Science and Technology. 
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